A&A manuscript no. 

(will be inserted by hand later ) 



Your thesaurus codes are: 

10.08.1;10.11.1;10.19.2;10.19.3;12.04.1;12.07.1 



ASTRONOMY 

AND 
ASTROPHYSICS 



Observation of Microlensing towards the Galactic Spiral 
Arms. EROS II 2 year survey * 



on • 
OS ■ 
t-h ; 

bJQ' 

3 : 

< 

o ■ 

m : 

m ■ 
> ■ 

On i 

o . 

(N ■ 

m ■ 

o : 

On . 
On ■ 

^ ' 

Or 

6 ■ 

a: 



F. Derue 1 , C. Afonso 2 , C. Alard 3 , J-N. Albert 1 , A. Amadon 2 , J. Andersen 4 , R. Ansari 1 , E. Aubourg 2 , 
P. Bareyre 5 , F. Bauer 2 , J-P. Beaulieu 6 , A. Bouquet 5 , S. Char 7 **, X. Chariot 2 , F. Couchot 1 , 
C. Coutures 2 , R. Ferlet 6 , J-F. Glicenstein 2 , B. Goldman 2 , A. Gould 8 ***, D. Graff 2 9 , M. Gros 2 , 
J. Ha'fssinski 1 , J-C. Hamilton 5 , D. Hardin 2 , J. de Kat 2 , A. Kim 5 , T. Lasserre 2 , E. Lesquoy 2 , C. Loup 6 , 
C. Magneville 2 , B. Mansoux 1 , J-B. Marquette 6 , E. Maurice 10 , A. Milsztajn 2 , M. Moniez 1 , N. Palanque- 
Delabrouille 2 , O. Perdereau 1 , L. Prevot 10 , N. Regnault 1 , J. Rich 2 , M. Spiro 2 , A. Vidal-Madjar 6 , 
L. Vigroux 2 , S. Zylberajch 2 
The EROS collaboration 
in memory of Sergio Char 



Laboratoire de l'Accelerateur Lineaire, IN2P3-CNRS, Universite de Paris-Sud, B.P. 34, 91898 Orsay Cedex, France 
CEA, DSM, DAPNIA, Centre d'Etudes de Saclay, 91191 Gif-sur-Yvette Cedex, France 
DASGAL, INSU-CNRS, 77 avenue de l'Observatoire, 75014 Paris, France 

Astronomical Observatory, Copenhagen University, Juliane Maries Vej 30, 2100 Copenhagen, Denmark 
College de France, LPCC, IN2P3-CNRS, 11 place Marcellin Berthelot, 75231 Paris Cedex, France 
Institut d'Astrophysique de Paris, INSU-CNRS, 98 bis Boulevard Arago, 75014 Paris, France 

7 Universidad de la Serena, Facultad de Ciencias, Departamento de Fisica, Casilla 554, La Serena, Chile 

8 Department of Astronomy, Ohio State University, Columbus, Ohio 43210, U.S.A. 

9 Department of Physics, Ohio State University, Columbus, Ohio 43210, U.S.A. 

10 Observatoire de Marseille, INSU-CNRS, 2 place Le Verrier, 13248 Marseille Cedex 04, France 

Received 13/03/1999, accepted 



Abstract. We present the analysis of the light curves 
of 8.5 million stars observed during two seasons by 
EROS (Experience de Recherche d'Objets Sombres), in 
the Galactic plane away from the bulge. Three stars have 
been found that exhibit luminosity variations compati- 
ble with gravitational microlensing effects due to unseen 
objects. The corresponding optical depth, averaged over 
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four directions, is r 



0.38; 



-0.53 
-0.15 



x 10 . All three can- 



1. Introduction 



didates have long Einstein radius crossing times (<~ 70 
to 100 days). For one of them, the lack of evidence for 
a parallax or a source size effect enabled us to constrain 
the lens-source configuration. Another candidate displays 
a modulation of the magnification, which is compatible 
with the lensing of a binary source. 

The interpretation of the optical depths inferred from 
these observations is hindered by the imperfect knowledge 
of the distance to the target stars. Our measurements are 
compatible with expectations from simple galactic models 
under reasonable assumptions on the target distances. 



Since the seminal paper of Bohdan Paczynski 1986, obser- 



vations have demonstrated that gravitational microlensing 
is an efficient tool to investigate the Milky Way struc- 
ture. After the first detections of microlensing effects to- 



wards the Large Magellanic Cloud (Alcock et al. 1993; 
Aubourg et al. 1993) and towards the Galactic bulge 



Send offprint requests to: M. Moniez: Moniez@lal.in2p3.fr; 
see also our WWW server at URL : 



(Udalski et al. 1993; Alcock et al. 1995), searches for mi- 
crolensing have entered an active era. The results of recent 
campaigns of observations are somewhat difficult to inter- 
pret. The negative search for short duration events and 
the rarity of long duration events found towards the Mag- 
ellanic Clouds (Ansari et al. 1996; Alcock et al. 1997a; 



Palanquc-Dclabrouille et al. 1998 and Alcock et al. 1998) 



http : //www. lal . in2p3 . f r/EROS/eros . html 
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imply that only heavy dark compact objects (M > 
10~ 2 M©) could account for a significant fraction (> 25%) 
of the halo mass required to explain the rotation curve 
of our galaxy. On the other hand, the optical depth mea- 
sured in the Galactic bulge direction ( Udalski et al. 1994 ; 
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Fig. 1. Expected optical depth (xlO 6 ) up to 7 kpc for the 
different components of the Milky Way as a function of 
the Galactic longitude at b — —2? 5 for two models of the 
Galaxy (see Table 4). The 4 directions towards the spiral 
arms are indicated. Notice that at this Galactic latitude 
the bulge and disc contributions are already reduced by 
a factor ~ 2 with respect to zero latitude; moreover, the 
bulge contribution also changes dramatically with the dis- 
tance to the target (here assumed to be at 7 kpc, while 
the Galactic Centre is at 8.5 kpc). 
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passbands. Photons are collected by two cameras, each 
equipped with a 2 x 4 mosaic of 2K x 2K LORAL CCDs 



Alcock et al. 1997b) appears to be significantly larger 
than expected from the contributions of the disc, the bulge 
and the halo. The latter measurements have led to the hy- 
pothesis of an ellipsoidal structure of the Galactic bulge. 
It is thus important to disentangle the different contribu- 
tions to this large optical depth. The estimate of the disc 
contribution can be refined by investigating lines of sight 
that do not go through the hypothetic ellipsoidal struc- 
ture (see Fig. [j}). Therefore the EROS team has chosen to 
search for microlensing not only towards the Magellanic 
Clouds and the Galactic bulge, but also in four regions 
of the Galactic plane, located at a large angle from the 
Galactic Centre. 

2. The observations 

2.1. Data taking 

Since July 1996 the EROS team has been using at La 
Silla observatory the MARLY telescope (1 m, f/5) with a 
dichroic beam-splitter allowing simultaneous imaging of a 
0.°7(cv) x l.°4(<5) field in EROS-visible and EROS-red wide 



(Bauer et al. 1997). In this analysis, seven sub-fields are 
considered per image because one of our 16 CCDs is not 
operational. The pixel size is 0.6 arcsec, and the median 
global seeing is 2 arcsec. CCDs are read-out in parallel 
in approximately 50 s, during which the telescope moves 
towards the next field to be imaged. After acquisition by 
the VME system, raw data are reduced by DEC-Alpha 
workstations using flat-field images taken at the begin- 
ning of the night. The DLT tapes produced are shipped to 
the CCPN (IN2P3 computing centre, CNRS) in Lyons, 
France, for subsequent processing. 

The two EROS passbands are nonstandard. EROS- 
red passband is centered on A = 762 nm with a full width 
half maximum AA ~ 200 nm, and EROS-visible passband 
is centered on A = 600 nm with AA ~ 200 nm. Our cal- 
ibration studies show that the corresponding magnitudes 
Reros and Veros are close to Cousins I and Johnson 
R to ±0.3 magnitudes. These Reros and Veros magni- 
tudes will be used all along this article. 



2.2. The targets 

Four different directions are monitored in the Galactic 
plane, away from the bulge, totalizing 29 fields which have 
a high stellar density, and cover a wide range of Galactic 
longitude. We refer to them as f3 & 7 Set, 7 Nor and 6 Mus. 
Table @ gives the coordinates and available data related 
to the monitored directions; Fig. || displays the positions 
of the fields in Galactic coordinates and Fig. || shows the 
observation periods and average time sampling. The ex- 
posure times of 2 and 3 minutes were chosen to optimize 
the global sensitivity of the photometric measurements 
taken during microlensing magnifications (a compromise 
between the number of measurements and their precision, 
see Mansoux 1997). 

By contrast with the Magellanic Clouds, the distance 
distribution of the monitored stars is imperfectly known, 
and should a priori vary with the limiting magnitude. In 
our detection conditions, the populations of stars used to 
obtain the optical depths given in Sect. 4.1 below are those 
described by the colour-magnitude diagrams of Fig. f§ An 
analysis of these diagrams has shown that their content 
is dominated by a population of source stars located ~ 7 
kpc away, undergo ing an interste llar extinction of about 
3 magnitudes (see Mansoux 1997 for more details). This 
distance estimate is in rough ag reement with the dista nce 
to the spiral arms d educed from Gcorgclin et al. 1994 and 
Russeil et al. 1998 and will be used in this paper. 



3. The search for lensed stars 

3.1. Data processing and analysis 

Light curves have been produced from the sequences of 
images using the specific software PEIDA (Photometrie 
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Fig. 2. Map of the Galactic plane fields (Galactic coordinates) monitored by EROS for the microlensing search. The 
shaded area represents the shape of the Galaxy. We have indicated our Galactic Centre fields and the 4 directions 
towards the arms. 
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Table 1. Description of the 29 fields monitored in the 
spiral arms program. This table gives the fields centres, 
the averaged number of measurements and the number of 
light curves analysed in this article. Field gn401 has not 
been studied yet. 



Field 


a (h:m:s) 


S (d:m:s) 


# of 


stars 


number 


J2000 


J2000 


meas. 


(million) 


Scutum (j3 Set) 


53 


1.96 


bs300 


18:43:22.0 


-07:40:53 


55 


0.33 


bs301 


18:43:27.0 


-06:13:42 


52 


0.23 


bs302 


18:46:16.0 


-07:22:45 


50 


0.32 


bs303 


18:46:20.0 


-05:55:35 


50 


0.28 


bs304 


18:49:21.0 


-06:45:51 


54 


0.43 


bs305 


18:52:26.0 


-06:35:44 


53 


0.37 


Scutum (7 Set) 


51 


1.70 


gs200 


18:28:03.0 


-14:51:06 


55 


0.32 


gs201 


18:31:15.0 


-14:14:38 


52 


0.30 


gs202 


18:31:33.0 


-12:48:53 


49 


0.29 


gs203 


18:34:22.0 


-14:31:39 


50 


0.34 


gs204 


18:34:28.0 


-13:04:31 


50 


0.45 


Norma (7 Nor) 


100 


3.01 


gn400 


16:09:45.0 


-53:07:03 


100 


0.37 


gn401 


16:18:22.0 


-51:44:43 






gn402 


16:14:57.0 


-53:04:35 


101 


0.25 


gn403 


16:22:28.0 


-52:06:20 


108 


0.30 


gn404 


16:19:09.0 


-53:26:38 


106 


0.29 


gn405 


16:26:52.0 


-52:21:02 


111 


0.22 


gn406 


16:23:54.0 


-53:43:53 


106 


0.33 


gn407 


16:31:31.0 


-52:28:44 


107 


0.22 


gn408 


16:28:42.0 


-53:51:58 


92 


0.22 


gn409 


16:15:51.0 


-54:48:45 


90 


0.26 


gn410 


16:20:30.0 


-55:04:18 


82 


0.23 


gn411 


16:09:37.0 


-55:10:07 


90 


0.32 


Musca (6 Mus) 


66 


1.77 


tm500 


13:27:04.0 


-63:02:18 


65 


0.31 


tm501 


13:31:18.0 


-63:34:41 


64 


0.33 


tm502 


13:34:52.0 


-64:10:30 


68 


0.33 


tm503 


13:23:58.0 


-64:59:52 


64 


0.29 


tm504 


13:12:12.0 


-64:06:49 


70 


0.23 


tm505 


13:16:15.0 


-64:40:50 


65 


0.28 


Total 


8.44 



et Etude dlmages Destinies a I'Astrophysique) , designed 
to extract photometric information in crowded fields 
( |Ansari 1996j) . Figure |^ shows the mean point-to-point rel- 
ative dispersion of the measured fluxes along those light 
curves as a function of Reros and Veros- A systematic 
search for microlensing events has been performed on a set 
of about 300,000 stars per field (8.5 x 10 6 light curves in 
total) , with an average of 70 measurements in the Veros 
and Reros colours (between July 1996 and November 
1997). 

This search is based on procedures and criteria similar 
to the ones deve loped for the EROS microlensing s earch in 
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the SMC data ( [Palanque-Delabrouille et al. 1998[ ). These 



criteria are designed to select microlensing events using 



20 30 40 50 60 70 80 90 20 30 40 50 60 70 80 90 100 

number of weeks since 01/01/96 

Fig. 3. Time sampling for each direction monitored to- 
wards the spiral arms, in number of measurements per 
week. 7 Nor is monitored between January and October, 
9 Mus is observed between January and August. (3 & 7 
Set are monitored between May and November. 



the expected characteristics of their light curves (one sin- 
gle peak, simultaneous in the visible and red bands, with 
a known shape), and to reject variable stars (rejection of 
stars lying in regions of the colour-magnitude diagrams 
mostly populated by variable stars, and requirement for 
the stability of the curve outside the peak). They are 
also tuned to be loose enough to avoid rejection of non- 
standard microlensing events, which have a different peak 
shape. 

The main difference from the SMC analysis is the re- 
jection against variable stars from the instability strip and 
against red giant variables. Unlike the Magellanic Clouds 
case where the positions of these populations in the colour- 
magnitude diagram are known a priori, the scatter and im- 
perfect knowledge of the distances and reddenings of our 
target stars make our colour-magnitude cut somewhat em- 
pirical. In particular, its acceptance is different from one 
field to another (see the excluded regions in Fig. |^) . 

3.2. The efficiency of the analysis 

To determine the efficiency of each selection criterion, we 
have applied them to Monte-Carlo generated light curves, 
obtained from a representative sample of the observed 
light curves, on which we superimpose randomly generated 
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Fig. 5. Relative frame to frame average dispersion of the 
luminosity measurements versus Re bos (upper panel) 
and Veros (lower panel), for stars with at least 50 re- 
liable measurements for each colour. This dispersion is 
taken as an estimator of the mean photometric precision. 
The superimposed hatched histograms show the magni- 
tude distribution of the stars in EROS bands. 



microlensing effects. The microlensing parameters are uni- 
formly drawn in the following intervals: impact parameter 
expressed in Einstein radius unit uq 6 [0,2], maximum 
magnification time in a research period T Q b s starting 50 
days before the first observations and ending 50 days after 
the last observations^], and Einstein radius crossing time 
At £ [1, 150] days. The analysis efficiency (or sampling 
efficiency) shown in Fig. ^ is relative to a set of unblended 
stars. Blending of a lensed star results in the mixing of 
unmagnified light with the expected light curve. The ef- 
ficiency of the search and the reconstruction of the lens- 
ing parameters are hence affected by this distortion. On 
the other hand, the set of stars that could undergo a de- 
tectable lensing effect is larger than the sample of the mon- 
itored stars, and the expected number of events is also af- 
fected. These effects, which are usually included as correc- 
tions to the sampling efficiency, are not taken into account 
in the present article and will require a specific study. For 
other targets (SMC, LMC), we have determined the cor- 
rection o f the detection rate due to blending to be les s than 
10-20% ( |Renault 1996] |Palanque-Dclabrouillc 19971) . 



1 Outside this domain of parameters, the detection efficiency 
is negligible. 




Event timescale (days) 

Fig. 6. EROS microlensing detection efficiency e(At) as a 
function of At, the Einstein radius crossing time. e(At) is 
the ratio of the number of simulated events, satisfying the 
selection criteria with duration At, with any uq and any 
date of peak magnification within the research period (650 
days), to the number of events generated with uq < 1. 



3.3. Results of the selection 

Three light curves satisfy all the requirements and are 
hereafter named candidates and labelled GSAl to 3. Fig- 
ures 0, H and [9] show the light curves of each candidate 
and Table § contains their characteristics. Measurements 
taken after Jan 1st, 1998 (date 2922) are shown, although 
they were not used in the selection. 



4. Optical depth and event timescales 

4-1. Optical depth 

The optical depth towards a pointlike source is defined 
as the fraction of time during which it undergoes a lens- 
ing magnification larger than 1.34. For a given target the 
measured optical depth r is computed from: 



1 



E 



At 



N obs T obs 2 ^ e(At) ' 

events x 7 



where N a b s is the number of monitored stars in the target, 
T b 8 is the duration of the search period (650 days for this 
2 year analysis) and e(At) is the average detection effi- 
ciency normalized to the microlensing events with impact 
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Table 2. Characteristics of the 3 microlensing candidates 



Candidate 


rjllUuz-LraAl 




Xt VJ O Z - O i\ O 


field 


7 Set 


7 Nor 


7 Nor 


Coordinates of star (J2000) 


a = 18h29m09.0s 


a = 16hllm50.2s 


a = 16hl6m26.7s 




S = -14° 15' 09" 


S = -52° 56' 49" 


S = -54°37'49" 


Galactic coordinates 


b = -If 67 


b= -If 17 


b = -2.°86 




I — 1 ( . ID 




/ QOO° QO 


Reros 


17.6 


17.7 


17.4 


Veros 


19.7 


19.4 


18.6 


Date of maximum magnification 


Aug. 3 rd , 1996 


Mar. 26 th , 1997 


Oct. 7 th , 1997 


Julian Day -21,447,891.5 


2408.6 ±0.1 


2642.9 ± 0.2 


2806.2 ± 1.1 


Einstein radius crossing 


73.5 ± 1.4 


98.3 ±0.9 


70.0 ±2.0 


time At = Re/Vt (days) 








Max. magnification 


26.5 ±0.6 


3.05 ± 0.02 


1.89 ±0.01 


Impact parameter (in Re) 


0.0378 ± 0.001 


0.342 ± 0.002 


0.593 ± 0.007 


X 2 of best fit 


185.7/163 d.o.f. 


551/425 d.o.f. 


445/427 d.o.f. 


Remarks 


M lens > 4.6 x 1O" 3 M 


binary source fit 






at 95% C.L. (see text) 


period 98 days (see text) 





Table 3. Contribution of the candidates to the optical 
depth t assuming the sources to be 7 kpc away. In the 
case of 9 Mus we give a 95% C.L. upper limit on the 
optical depth contribution from events with At = 80 days 
(e(80 days) = 18%). 



Target 


Scutum 


Norma 


Musca 


Direction 


P Set 7 Set 


7 Nor 


6 Mus 


<6> 


-2f5 -2.°6 


-2.° 7 


-l.°8 


< I > 


27.°0 18.°6 


331.° 2 


306.° 4 


Events detected 


none GSA1 


GSA2 & 3 


none 


At (days) 


73 


98 & 70 


(80) 


e(Ai)(%) 


10 


27 


(18) 


r(event)(xl0 6 ) 


1.02 


0.29 & 0.21 




f (target) (xlO 6 ) 


0.47 


0.5 


< 1.82 at 








95% CL 


t averaged over the 4 directions (xl0°) 


0.38^.i5 



parameter uq < 1, whose maximum magnification takes 
place within the research period. The contribution of the 
candidates to the optical depth is given in Table 3. 

We have modeled the Galaxy in two different ways 
using three components: a central bulge, a disc and a dark 
halo. The density distributio n for the bulge - a barlike 
triaxial model - is taken from Dwek et al. 1995 model G2 
(in cartesian coordinates): 



M 



PBulge 



Bulge e ~r 2 /2 



8irabc 



where Msuige is the bulge mass, and a, b, c are the length 
scale factors. The bar major axis is inclined at an angle of 
15° with respect to the Sun-Galactic Centre line. 



We compare our measurements with the predictions 
of two models with extreme disk contributions (here the 
main structure involved in microlensing) . The first one has 
a "thin" disc and a standard isotropic and isothermal halo 
(Model 1) with a density distribution given in spherical 
coordinates by: 



P(r) = Pho 



Rq + R c 
r 2 + R?. 



where phQ is the local halo density, Rq — 8.5 kpc is the 
distance between the Sun and the Galactic Centre, and 
R c = 5 kpc is the Halo "core radius". The matter dis- 
tribution in the disc is modeled in cylindrical coordinates 



by a double exponential (see e.g Bienayme et al. 1987 and 
|Schaeffer et al. 1998| ): 



Pthin(R, z) 



-'thin 



2H, 



thi 



■ exp 



<R-Re) 

Rthin 



exp 



1*1 



H, 



l hi 



Tithin is the column density of the disc at the Sun position, 
Hthin is the height scale and Rthin is the length scale of 
the disc. 

The second model (Model 2) has a "thin" and a "thick" 
disc, and a very light halo. Both models share the same 
bulge contribution. The model parameters are summa- 
rized in Table 4. Fig. [l] shows the expected optical depth 
up to 7 kpc as a function of longitude for both models, at 
the average latitude of our fields b = —2.5°. As the main 
contribution comes from the thin disc (about 90%), vari- 
ations of the optical depth from field to field due to the 
range of 2 to 3° in latitude can reach ~ 30% in the case 
of 7 Nor, and ~ 20% for the other targets. 

The expected optical depth, averaged over the four di- 
rections, is 0.60 x 10~ 6 for model 1, and 0.70 x 10~ 6 for 



F. Derue et al. : Observation of Microlensing towards the Galactic Spiral Arms 



7 









EROS2- 


GSA1 (R E ros) 








25 




ft 




- 




20 
1 5 


: | 


\ 








10 

5 


/ 


\ 

V 












_ 




-j 


i 


O — ' — ' — ' — ' — >- 
23SO 2400 


2420 2440 2460 24SO 2500 




\ 

M 






. --- u - 








f • 








2400 2500 2600 2700 2800 2900 3000 31 OO 3200 

Time (JD — 244789 1 .5) 



2400 2500 2600 2700 2800 2900 3000 31 OO 3200 

Time (JD — 244789 1 .5) 





I 




EROS2 — GSA1 (V EROS ) 




j 


25 


A 










20 












1 5 


j 

/ 










10 

5 


/ 
/ 


\ 

V 






I 










- 

- 1 


t 


23SO 2400 


2420 2440 2460 24SO 2500 




* 








-# 






















0.5 



2400 2500 2600 2700 2800 2900 3000 31 OO 3200 

Time (JD — 244789 1 -5) 



2400 2500 2600 2700 2800 2900 3000 31 OO 3200 

Time (JD — 244789 1 .5) 



Fig. 7. Magnification curves of the microlensing candidate 
GSA1 in the direction of 7 Set. The fitted standard mi- 
crolensing curve is superimposed (solid line). 



model 2. These estimates vary by 50% if the average dis- 
tance of the sources is changed by 2kpc, or if the parame- 
ter Yjthin is changed by 25M /pc 2 . The measured optical 
depth, averaged over the four directions is: 

f = 0.38±g:fi x 1CT 6 . 

The confidence interval reported here takes into account 
Poisson fluctuations and the possible event timescale vari- 
ations inside the range [71,98] days. The comparison of 
Fig. [I] with Table 3 also shows that the measured opti- 
cal depths in the four directions are compatible with the 
predictions of both models. 

4-2. Microlensing event timescales 

The duration of the three events is long (~ 80 days in 
average). Fig. |l^ shows the expected event duration dis- 
tribution towards 7 Nor within the framework of model 1. 
This distribution is obtained assuming the following mass 
functions and kincmatical characteristics: 



Fig. 8. Magnification curves of the microlensing candidate 
GSA2 in the direction of 7 Nor. The solid line shows the 
fitted microlensing curve taking into account the modula- 
tion due to a dominant source orbiting in a binary system 
with period Pq = 98 days. The dashed line corresponds to 
the best standard microlensing fit. 



— We assume that lenses belonging to the (non-rotating) 
halo have the same mass (0.5M Q ); their velocities 
transverse to the line of sight of disc stars fol- 
low a Boltzmann distribution with a dispersion of 
~ 150km/s. Thus the expected duration of microlens- 
ing events is small (see Fig. |To| ) . 

— For lenses belonging to the bulge, the mass function is 
taken from Richer et al. 1992 and the velocities trans- 



verse to the line of sight of disc stars also follow a Boltz- 
mann distribution with a dispersion of ~ HOkm/s. The 
expected rate of microlensing events due to this struc- 
ture towards the direction considered in Fig. |l^ ('-/Nor) 
is found to be negligible. 

The disc lenses mass function is taken from 



Gould 1997, which is derived from HST observations. 



Disc lenses are subject to a similar global rotatio n as 
the observer and the sources ( Brand fc Blitz 1993 ). We 
assume a negligible particular motion of the monitored 
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Fig. 9. Magnification curves of the microlensing candidate 
GSA3, in the direction of 7 Nor. 



sources with respect to the spiral ar ms, because they 

the 



are probably young stars. Following |Gricst 1991 , 
motion of the Sun relative to the Local Standard of 
Rest is taken as (i>©_r = 9, v Q g = 11, v Qz = 16) (in 
km/s) ; the velocity dispersions of the lens population 
are expected to be (a(Vn) = 40, a(Vg) = 30, <r(V z ) — 
20) (in km/s). As the disc lenses have a low velocity 
relative to the line of sight, disc-disc events have longer 
timescales, as can be seen on Fig. [l(| 

Assuming that all three lenses belong to the halo leads 
to large probable masses (> 2M Q ). As such high mass 
lenses would be visible stars that we do not observe (dis- 
missing the unlikely possibility that the three could be 
neutron stars), the disc-disc lensing hypothesis is more 
probable. This hypothesis could also explain the fact that 
4 events over the 45 found by the MACHO Collabora- 



tion towards the Galactic Centre (see Alcock et al. 1997b) 
have a timescale > 50 days, significantly larger than the 
mean for the whole sample (21 days). Indeed, the contri- 
bution of these 4 events to the total optical depth is about 
0.6 x 10 -6 , which is compatible with the expectation of 
the disc lensing contribution (~ 0.7 x 10~ 6 for Model 1). 



Table 4. Parameters of the galactic models used in this 
article. 





Parameter 


Model 1 


Model 2 


Bulge 


(kpc) 
b (kpc) 
c (kpc) 

M bu i ge (xW 10 M Q ) 


1.49 
0.58 
0.40 
2.1 


Thin disc 


E (M pc-*) 
H (kpc) 
R (kpc) 

Mthin(XW Mq) 


50 
0.325 
3.5 
4.3 


Thick disc 


E (M Q pc-^) 
H (kpc) 
R (kpc) 

M tMcfc (xlO lo M ) 




35 
1.0 
3.5 
3.1 


Halo 


p h& (M pc- J ) 
Rc (kpc) 

M in 60 kpc (1O 1O M ) 


0.008 
5.0 
51 


0.003 
5.0 
7 


Predictions 


p Q (M pc~ J ) 
Vrot at sun (km s _1 ) 
Vrot at 20 kpc 
Vrot at 60 kpc 


0.085 
211 
203 
200 


0.098 
222 
180 
140 



5. Detailed analyses of the candidates 

Detailed analyses have been performed on GSAl and 
GSA2 candidates, whose light curves are sufficiently sam- 
pled. For this purpose, we have used all available data on 
the candidates (i.e. three years of observations). 



5.1. Candidate EROS2-GSA1 

Candidate GSAl exhibits a large magnification (A pea k > 
25 at 95% C.L.) with no detectable b lendingp| and no de- 
tectable parallax effect ( |Gould 1992 ). Despite the small 
impact parameter, no evidence for a distortion of the curve 
due to the non-zero size of the lensed source is found. A 
fitting procedure allows one to put an upper limit to the 
ratio of the angular stellar radius 9+ to the angular Ein- 
stein radius 9e of the lens: 



< 0.066 at 95% C.L 



4GM 1 
—pi 



xD 



where M is the lens mass, D the distance of the observer 
to the source and xD its distance to the lens. From the 
position of the source in the colour-magnitude diagram we 
can assume that its temperature is comparable to or lower 
than that of the Sun. We then obtain: 



00 at 10 pc 



> 



10< 



EROS0 



-v° 



2 Fits including more than 10% blending are excluded at 95% 
C.L. Moreover, since the magnification curves are the same in 
both colours, the amplified star and an hypothetical blending 
star would have to have -by chance- the same colour. 
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This implies a lower limit on the projected speed of the 
lens onto the local transverse plane: 

v = r E jAt > 31.6 km/s at 95% C.L. 

Since typical values expected for this quantity are of the 
order of v ~ 100 km/s, it is clear that this is not a very 
constraining limit either. Fig. ^ shows the excluded area 
for the tip of the 2-dim. vector v in the local trans- 
verse plane. The lower limit we find for the modulus v 
(and te) is relatively small, because the best parallax 
fit, which is not significantly better than a standard fit 
(183.7/161 d.o.f. compared with 185.7/163 d.o.f.), is ob- 
tained for v = (—17.1, —31.1) km/s in the frame of Fig. 
|TT| . This special configuration produces a distorted light 
curve which diverges only marginally from the standard 
fitted curve, and only during periods where the measure- 
ments are not very precise. 

From the definitions of 9e and ?e one gets the relation: 



"0 20 40 60 80 100 120 140 

Event timescale At (days) 

Fig. 10. At distribution of the events expected during the 
search period T Q b 8 per 10 7 monitored stars towards 7 Nor, 
within model 1. This distribution takes into account the 
detection efficiency e(At). The durations of the three can- 
didates are marked. Notice that the bulge contribution is 
negligible. 



where the apparent magnitude V^ os ^ should be cor- 
rected for the interstellar extinction. Ignoring this correc- 
tion leads to the conservative limit t?* > 2.37 x 10~ 12 rad, 



implying that E > 3.6 x 10" 11 
proper motion of the deflector: 



rad and that the angular 



u (km/s/kpc) = 3.57 x 10 11 H^f^K > 0.17 at 95% C.L. 
r v ' ' ; At (days) 

Assuming the interstellar extinction to be 3 magnitudes 
(see Sect. 2.2) leads to > 9.42 x 10~ 12 rad, implying 
that 9 E > 1.43 x 10~ 10 rad and: 

H > 0.69 km/s/kpc at 95% C.L. 

Since the typical expected value (based on Galactic kine- 
matics) is [i ~ 4 km/s/kpc, this limit does not probe the 
expected range of parameter space. 

On the other hand, the standard microlensing fit is 
satisfactory, and taking into account the possibility of a 
parallax effect does not result in a significantly better fit. 
We can then infer a lower limit on f E , the size of the 
Einstein ring projected from the source onto the observer's 
plane: 



—vfiAt 2 



9 E xD 



> 1.33 AU at 95%CL. 



from which we derive a lower limit on the lens mass by 
combining the two constraints from finite size and parallax 
analysis, i.e. at 95% C.L.: 

M > 1.2 x 1O~ 3 M ignoring interstellar extinction, 
M > 4.6 x 10 _3 M Q with 3 magnitudes of extinction. 

Note that this limit is independent of the lens and source 
distances to the observer. The fact that this lower limit is 
so much smaller than the expected mass for such lenses 
follows immediately from the fact that each limit from 
which it is derived is not probing the expected range of 
kinematic parameters. Fig. [T^ shows the excluded areas 
in the M versus x plane, from the finite size study and 
the parallax analysis, assuming the source to be located 7 
kpc away. Exclusion curves for the two hypotheses on the 
interstellar absorption are shown. A better characteriza- 
tion of the lensed source should allow one to refine these 
preliminary studies. 

Finally, a limit on the lens luminosity can be derived 
from the maximum blending limit allowed by the fit. At 
95% C.L., the apparent magnitude of the lens is at least 
2.5 magnitudes above the measured baseline magnitude, 
i.e. Reros > 20.1 and V EROS > 22.2. 

5.2. Candidate EROS2-GSA2 

GSA2 has residuals to a standard microlensing fit which 
clearly exhibit a modulation; a period of ~ 54 days is 
found in the residuals of the standard fit during the magni- 
fication. On the other hand, taking into account a possible 
parallax effect does not significantly improve this fit. Fur- 
thermore, we know that the lensing of a periodic variable 
star would produce a modulation of the light curve whose 
amplitude should follow the magnification. As we do not 
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Fig. 11. Excluded area for the tip of the 2-dim. vector 
v in the local transverse plane. The smallest speed com- 
patible with the observations corresponds to a deflector 
with a velocity oriented towards the lower spike. The ar- 
row shows the earth projected velocity at the maximum 
magnification time. The straight line is the intersection of 
the Galactic plane. The projection of the earth trajectory 
is indicated by the open circles, starting 60 days before 
maximum, ending 70 days after maximum, with 10 days 
spacing (right scale in AU). NEP is the North Ecliptic 
Pole direction, NGP is the North Galactic Pole direction. 



detect a significant modulation in the non-magnified part 
of the light curve, the most probable origin of this modula- 
tion is that one dominant luminous source or two sources 
orbit around the centre of gravity of a binary system, in- 
ducing a wobbling of the line of sight with respect to the 
trajectory of the lens 0. In this configuration, a modula- 
tion of the light curve is expected to occur only during the 



Thi s type of configuration has been s tud- 



ied by |Gricst fc Hu 1992| , |Sazhin fc Cherepashchuk 1994| and 
Han & Gould 1997, and w as already mentioned in our earlier 
article (Ansari et al. 1995). 



10 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x = D rt /D. 



Fig. 12. 95% C.L. exclusion diagram from the parallax 
analysis and from the finite source size study for the can- 
didate EROS2-GSA1. The lower curve in the right part 
is the conservative limit ignoring interstellar absorption, 
the higher curve corresponds to a 3 magnitude absorp- 
tion. These curves assume the source to be located 7 kpc 
away. If the source is assumed to be at some other distance 
D, one has to scale the increasing functions by the factor 
D/(7 kpc), and the decreasing function by (7 kpc)/£>. 



lensing peak. We obtain equally satisfactory fits for two 
classes of circular binary systems (with at best x 2 = 551 
for 425 d.o.f., to be compared with 876 for 439 d.o.f. for a 
standard microlensing fit). 

The first class includes systems with a dominant source 
orbiting with a period P ~ 51 days, and with the pro- 
jected semi-axis p — ax/ Re ~ 0.04 The second class 
corresponds to systems with two luminous stars orbiting 
with a period P ~ 98 days, with luminosity and mass 
ratios around 1/3-1/2, and a projected distance between 
the two components of about p = ax /Re — 0.4. Given 
the domain of values for x and Re, these parameters are 
compatible with physically acceptable masses of the bi- 
nary satisfying Kepler's third law. Many more configura- 
tions can fit the observations if we consider systems with 
elliptic orbits. A spectroscopic study of the source is under 
way in order to test its binarity. If the source proves to be 



4 Fits which are slightly less good can be obtained with 
larger periods, but in this case only the first shoulder (around 
date=2600) is correctly matched. 
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Fig. 13. Relation between Mi ens x D and x assuming a 
binary system with P = 100 days, p — 0.4 and a total 
mass M T of 1M© or 10M©. 



a spectroscopic binary, then one can get a good estimate 
or constraint on the angular Einstein radius 9e (and then 
on the angular proper motion of the lens p). 

Expressing the Einstein radius and Kepler's third law 
in the relation p — ax/ Re leads to: 



Mi 



lens 



D 



M 7 kpc 



6.7.10' 



6 


Po 


4 
3 


'M T ~ 




1 day_ 




[M e \ 



1 



Fig. [13] illustrates this relation between M\ ens x D and x 
assuming the source to be a binary system with P — 100 
days, with p = 0.4 and with two different hypotheses for 
the mass Mt of the system. 

6. Conclusion 

We have searched for microlensing events with durations 
ranging from a few days to a few months in four Galac- 
tic disc fields lying 18° to 55° from the Galactic Cen- 
tre. We find three events that can be interpreted as mi- 
crolensing effects due to massive compact objects. Their 
long duration favours the interpretation of lensing by ob- 
jects belonging to the disc instead of the halo. The aver- 
age optical depth measured towards the four directions is 
f = 0.38to^5 x 10~ 6 . Assuming the sources to be 7 kpc 
away, the expected optical depths from two different galac- 
tic models vary from 0.60 to 0.70 x 10~ 6 , in agreement with 
our measurement. 



One event displays a modulation of the magnification 
which is compatible with the lensing of a binary source. 
More information about the configuration of this possi- 
ble binary source, and about the lens mass acting on the 
strongly amplified source, are expected from complemen- 
tary observations. No evidence for parallax and blending 
effects has been found. 

The observations continues towards the spiral arms. 
More accurate measurements should be obtained with the 



increase of statistics ( Derue 199E ), allowing one to esti- 
mate the disc contribution to the optical depth towards 
the bulge and the Magellanic Clouds. 
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